Testing the Chiral Magnetic Effect with Central U + U colhsions 
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A quark interaction with topologically nontrivial gluonic fields, instantons and sphalerons, violates 
V and CV symmetry. In the strong magnetic field of a noncentral nuclear collision such interactions 
lead to the charge separation along the magnetic field, the so-called chiral magnetic effect (CME). 
Recent results from the STAR collaboration on charge dependent correlations are consistent with 
theoretical expectations for CME but may have contributions from other effects, which prevents 
definitive interpretation of the data. Here I propose to use central body-body U + U colhsions 
to disentangle correlations due to CME from possible background correlations due to elliptic flow. 
Further more quantitative studies can be performed with collision of isobaric beams. 
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I. INTRODUCTION 

Quantum cliromodynamics (QCD), the theory of 
strong interactions, is a non-Abelian gauge theory that 
possesses multiple vacua characterized by Chern-Simons 
numbers. QCD links chiral symmetry breaking and the 
origin of hadron masses to the existence of topologi- 
cally nontrivial classical gluonic fields, instantons and 
sphalerons, describing the transitions between the vac- 
uum states with different Chcrn-Simons numbers. Quark 
interactions with such fields change the quark chirality 
and are V and CV odd. For a review, see Refs. [U [2]. 
Though theorists have little doubt in the existence of 
such fields, they have never been observed directly, e.g. 
at the level of quarks in the deep inelastic scattering. It 
was suggested in Ref. [3] to look for metastable V and 
CV odd domains, space-time regions occupied by a clas- 
sical field with a nonzero topological charge, in ultrarel- 
ativistic heavy ion collision. Some earlier discussion of 
how one could observe this local strong parity violation 
can be found in Refs. [3H5]. 

The situation with experimental search for the local 
strong parity violation drastically changed once it was 
noticed [51 [7] that in noncentral nuclear collisions it would 
lead to the asymmetry in the emission of positively and 
negatively charged particle perpendicular to the reaction 
plane. Such charge separation is a consequence of the dif- 
ference in the number of quarks with positive and neg- 
ative helicities positioned in the strong magnetic field 
(~ 10"'^^ T) of a noncentral nuclear collision, the so- 
called chiral magnetic effect (CME) [6l |8]. The same 
phenomenon can also be understood as an effect of the 
induced electric field that is parallel to the static external 
magnetic field, chiral magnetic induction, which occurs in 
the presence of topologically nontrivial gluonic fields [5]. 
It has been also argued that the charge separation could 
have origin in nonzero vorticity of the system created in 
noncentral collisions fT . Chiral magnetic effect has been 
observed in the lattice QCD calculations [TDHT^ . Newer 
developments in this field has been recently discussed at 
the RIKEN BNL workshop [H]. 

An experimental observation of CME would be a direct 



proof for the existence of topologically nontrivial vacuum 
structure and would provide an opportunity for a direct 
experimental study of the relevant physics. The diffi- 
culty in experimental observation of CME comes from 
the fact that the direction of the charge separation varies 
in accord with the sign of the topological charge of the 
domain. Then the observation of the effect is possible 
only by correlation techniques. According to Refs. [HHH] 
the charge separation could lead to asymmetry in parti- 
cle production (7V_ - 7V+)/(iV_ + N+) ~ Q/N„+, where 
Q — 0, ±1,±2,... is the topological charge and is 
the positive pion multiplicity in one unit of rapidity - 
the typical scale of such correlations. It results in corre- 
lations of the order of 10~*, which is accessible in current 
high statistics heavy ion experiments. An observable di- 
rectly sensitive to the charge separation effect has been 
proposed in Ref. [14] . It is discussed in more detail below. 

Recent STAR results [HI [16] on charge separation rela- 
tive to the reaction plane consistent with the expectation 
for CME can be considered as evidence of the local strong 
parity violation. The ambiguity in the interpretation of 
experimental results comes from possible contribution of 
(the reaction plane dependent) correlations not related to 
CME. As the detailed quantitative predictions for CME 
do not yet exist, it is difficult to disentangle different 
contributions. A key ingredient to CME is the strong 
magnetic field, while the background effects originate in 
the elliptic flow. In noncentral collisions of spherical nu- 
clei such as gold, both, magnetic field, and the elliptic 
flow are strongly (cor)related to each other. In order to 
disentangle the two effects one has to flnd a possibility to 
signiflcantly change the relative strengths of the magnetic 
fleld and the elliptic flow. The discussion of such a pos- 
sibility provided hy U + U collisions is the subject of this 
Letter. Necessary quantitative estimates are obtained 
using Glauber Monte-Carlo simulations. I estimate the 
magnetic field following the approach of Ref. [17 . Esti- 
mates of elliptic flow are based on the assumption that it 
scales with initial (participant) eccentricity of the nuclei 
overlap region. 

In noncentral nuclear collisions particle distribution in 
azimuthal angle is not uniform. The deviation from a 
flat distribution is called anisotropic flow and often is 
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described by the Fourier decomposition [151 US] (for a 
review, see Ref. [20]): 

— oc 1 + 2wi Q cos(A(?!)) + 2i;2 Q cos(2A(/)) + ... 
dq) 

+ 2ai,Qsin(A0) + 2a2,Qsin(2A(/)) + ... , (1) 

where A0 — {(j) — rp) is the particle azimuth relative 
to the reaction plane, vi and V2 account for directed and 
elliptic flow. Subscript a denotes the particle type. Be- 
cause of the "up-down" symmetry of the collisions a„ co- 
efficients are usually omitted. CME violates such a sym- 
metry. Although the "direction" of the violation fluctu- 
ates event to event and on average is zero, in events with 
a particular sign of the topological charge, the average is 
not zero. As a result, it leads to a nonzero contribution to 
correlations. One expects that the first harmonic would 
account for the most of the effect. To measure (ai.Qai,^), 
it was proposed [TJ] to use the correlator: 

(cos(0„ + 0/3-2«'fl,p)) = (2) 
(cosA(/)q cos A0^) — (sinA(/)a sinA^^j) 

= [{vi.aVl.ll) + Bin] - [{aq^aO-l^p) + Bout] 

^ — {ai,aai, is) + [Bin — Bout]- (3) 

This correlator represents the difference between corre- 
lations "projected" onto the reaction plane and the cor- 
relations projected onto an axis perpendicular to the re- 
action plane (a more detailed discussion of that can be 
found in Refs. [HKTB]. The key advantage of using such 
a difference is that it removes the correlations among 
particles a and /3 that are not related to the reaction 
plane orientation. The remaining background in Eq. [3] 
Bin — Bout, is due to processes in which particles a and (3 
are products of a cluster (e.g. resonance, jet, di-jets) de- 
cay, and the cluster itself exhibits elliptic flow or decays 
(fragments) differently when emitted in-plane compared 
to out-of-plane. The corresponding contribution to the 
correlator can be estimated as 1 14, 16 : 



(cos(0„ + 0^-2«'fl,p)) 
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where {■■■)ciust indicates that the average is performed 
only over pairs consisting of two daughters from the same 
cluster. This kind of background can not be easily elimi- 
nated or suppressed and constitute the main uncertainty 
in the interpretation of the STAR results. To address its 
contribution one has to rely on model calculations [16 . A 
better approach would be to perform experiments where 
the relative contribution of CME and background can 
be varied. Note that the background is proportional to 
the elliptic flow present in the event, V2^ciust in Eq, 
It is not clear how one could suppress elliptic flow an* 
at the same time preserve strong magnetic field needed 
for CME. But the opposite, collisions with strong ellip- 
tic flow and no (or almost no) magnetic field, seems to 
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FIG. 1: Schematic view of central U + U collisions: (a) tip-tip 
and (b) body-body. 



be possible. This can be achieved in central body-body 
U + U collisions. Uranium nuclei are not spherical and 
have roughly ellipsoidal shape. Central collision, when 
most of the nucleons interact, can have different geome- 
try, ranging from the so-called tip-tip collisions to body- 
body collisions [2T], see Fig. [l] Unlike tip-tip collisions, 
body-body ones would exhibit strong elliptic flow. Nei- 
ther would lead to a strong magnetic field; consequently, 
a very weak signal due to CME is expected, while back- 
ground would be much stronger in body-body configura- 
tion compared to tip-tip configuration. 

Collisions of uranium nuclei were first proposed for 
RHIC by P. Braun-Munzinger [22] with the goal to 
achieve higher energy density compared to Au-|-Au colli- 
sions. The idea was later elaborated in Refs. [231 [H], in 
particular pointing out an important possibility to study 
elliptic flow at such high energy densities. 

At RHIC one can select central collisions by requir- 
ing low signal in the zero degree calorimeters that detect 
spectator neutrons. Below I discuss how one might "con- 
trol" the geometry of the collision, and, consequently, the 
relative strengths of the signal due to CME and back- 
ground, by selecting events based on multiplicity, signal 
in the zero degree calorimeters, and the magnitude of the 
flow vectors. 

The magnetic field strength at a position r and time t 
is defined by the Lienard-Wiechert potentials 



eB{t,r) 



ttEM 



(1 



(Rn — RnVn) 



X R,, 



(5) 



where ckem ~ 1/137 is the fine-structure constant, and 
e„ is the electric charge of the nth particle in units of 
the electron charge. R„ = r — r„, where r„ is the radius 
vector of particle, v„ is particle velocity. The quantities 
v„ and r„ are taken at retarded time t' = t — |r — r„(t')|. 
Summation runs over all charged spectators. Spectator 
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contribution to the magnetic field is dominant at early 
times [3 [17] ; we also use this approximation in our esti- 
mates. Because of the Lorentz contraction, in collisions 
of ultrarelativistic nuclei, the longitudinal size of the nu- 
cleus is negligible compared to the transverse size, and 
the time dependence of the magnetic field is totally deter- 
mined by the gamma-factor (energy) of colliding nuclei. 
We are interested only in a relative change in the strength 
of the magnetic field in collisions at different centralities 
and different configurations. For this, it is sufficient to 
calculate the magnetic field only at t = (the time the 
two nuclei collide). At ^Jsnn = 200 GeV, the collision 
energy used in our estimates, the magnetic field aX t = 
is about factor of 2 smaller compared to the maximum 
value reached approximately at i « 0.05 fm/c |17) . 

Elliptic flow is determined by the geometry of the over- 
lapping zone. We assume V2 — K£p, where Sp is the so- 
called participant eccentricity. We consider only events 
within < 5% of the most central collisions, for which 
K ~ const. For definitions of eccentricity and details 
of experimental measurements of V2/S, see review |20j . 
We take k ~ 0.2 based on experimental measurements. 
Initial eccentricity, magnetic field, and charged particle 
multiplicity at midrapidity are calculated using Glauber 
Monte-Carlo with all parameters taken the same as used 
in Ref. 

The effect of nonsphericity of uranium nuclei is clearly 
visible in Fig. [2] which shows the distribution of events 
in elliptic flow V2 in event samples with number of spec- 
tators, Nsp < 20. The average elliptic flow is almost a 
factor of 2 larger in U + U collisions compared to Au-t-Au, 
which would mean a strong increase in the background 
correlations compared to that of due to CME. The re- 
quirement of the same number of spectators assures that 
the magnetic field is not very different in the two systems; 
it is slightly lower in U + U collisions [see Fig. [3]d]. The 
condition Ngp < 20 selects about 1.5% of the most cen- 
tral events in U + U collisions and about 2.3% in Au-|-Au 
collisions. 

Within the event sample selected on the basis of num- 
ber of spectators, it would be instructive to study the 
dependence of the signal on the magnitude of the elliptic 
flow. As a measure of the latter we use the magnitude of 
the flow vector q = Q2/VM, where 
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Q2,x = ^ cos(20i), Q2^y = ^ sin(2(/)i). 



(6) 



and the sum runs over all particles in a given momentum 
window. We calculate the flow vector based on charged 
particle multiplicity in 2 units of rapidity. As shown in 
Fig. [3) the elliptic flow is strongly correlated with g, and 
at the same time the magnetic field has almost no q- 
dependence. It means that the correlator, Eq. |3j used to 
measure the signal would stay constant if the signal is 
mostly determined by CME, and increase strongly with 
q if it is due to the background effects. For such a test 
U + U collisions provide a significantly better opportu- 
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FIG. 2: (Color online) Event distributions in 112 for Au+Au 
and U + U collisions in event samples with the number of 
spectators Nsp < 20. 
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FIG. 3: (Color online) Elliptic flow and the magnetic field (in 
arbitrary units) as a function of q, the magnitude of the flow 
vector, in events with the number of spectators Nsp < 20. 
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FIG. 4: (Color online) Elliptic flow and the magnetic field 
(arbitrary units) in Au+Au and U + U collisions as function 
of multiplicity. The arrows indicate the multiplicities corre- 
sponding to the top 2% of the collision cross section. 

nity than Au+Au collisions. First, the relative variation 
in V2 is almost a factor of 2 larger that that in Au+Au 
collisions. Also, the variation in elliptic flow in Au+Au 
collision is mostly determined by fluctuations in the ini- 
tial eccentricity, which are still not very well known. In 
U+U collisions the elliptic flow variation is mostly due to 
variation in orientation of the nuclei at the moment of col- 
lision. The corresponding estimates have much smaller 
uncertainty. 

While selection of the events based on the number of 
spectators is very useful, it seems to be also possible to 
disentangle CME and background correlations based only 
on the dependence of the signal on charged multiplicity. 
Figure |4] presents the dependence of the elliptic flow and 
magnetic field on charged multiplicity. The elliptic flow 
dependence is different for two systems, with U + U col- 



lisions exhibiting a characteristic kink (cusp) at multi- 
plicity ^ 1000 ^25J, reflecting the fact that high(er) mul- 
tiplicity events have predominantly tip-tip orientation; 
the latter also leads to a decrease in elliptic flow. Being 
mostly determined by correlation of the multiplicity with 
the number if participants, the magnetic field has simi- 
lar dependence on multiplicity for both collision systems. 
The difference in the dependencies of the magnetic field 
and elliptic flow on charged multiplicity can be used a as 
a test for the nature of correlations contributing to the 
signal. 

The charge separation dependence on the strength of 
the magnetic field can be further studied with collision 
of isobaric nuclei, such as and ^oZr. These nuclei 

have the same mass number, but differ by the charge. 
The multiparticle production in the midrapidity region 
would be affected very little in collision of such nuclei, 
and one would expect very similar elliptic flow. At the 
same time the magnetic held would be proportional to 
the nuclei charge and can vary by more than 10%, which 
can results in 20% variation in the signal. Such variations 
should be readily measurable. The collisions of ^^^Ru and 
40 isotopes have been successfully used at GSI [55] in a 
study of baryon stopping. Collisions of isobaric nuclei at 
RHIC will be also extremely valuable for understanding 
the initial conditions, and in particular the initial velocity 
fields, the origin of directed flow, etc. 

In summary, the estimates presented in this Letter 
show that a detailed analysis of central Au+Au and U+U 
collisions should be able to disentangle CME and back- 
ground correlations contributing to the signal observed 
by STAR. 
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